NACA TN 2931 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2931 


A METHOD FOR DETERMINING CLOUD-DROPLET IMPINGEMENT 

ON SWEPT WINGS 

By Robert G. Dorsch and Rinaldo J. Brun 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


NACA 


Washington 

April 1953 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2931 


A METHOD FOR DETERMINING CLOUD-DROPLET IMPINGEMENT ON SWEPT WINGS 
By Robert G. Dorsch and Rinaldo J. Brun 


SUMMARY 

The general effect of wing sweep on cloud-droplet trajectories 
about swept wings of high aspect ratio moving at subsonic speeds is 
discussed. A method of computing droplet trajectories about yawed 
cylinders and swept wings is presented, and illustrative droplet tra- 
jectories are computed. 

A method of extending two-dimensional calculations of droplet 
impingement on nonswept wings to swept wings is presented. It is shown 
that the extent of impingement of cloud droplets on an airfoil surface, 
the total rate of collection of water, and the local rate of impinge- 
ment per unit area of airfoil surface can be found for a swept wing 
from two-dimensional data for a nonswept wing. The impingement on a 
swept wing is obtained from impingement data for a nonswept airfoil 
section which is the same as the section in the normal plane of the 
swept wing by calculating all dimensionless parameters with respect to 
flow conditions in the normal plane of the swept wing. 


INTRODUCTION 

The efficient design of a thermal ice-protection system for a wing 
and the assessment of the heat requirements of the system necessitate 
a knowledge of the expected distribution of impinging supercooled cloud 
droplets over an airfoil surface exposed to a given set of meteorological 
conditions in flight. Because a closed-form solution for the equations 
of motion of a droplet in the flow field around an airfoil generally 
cannot be obtained, the impingement calculations are best performed with 
the aid of a differential analyzer or other computing machine. Two- 
dimensional calculations of droplet impingement on various airfoils, 
obtained with computing machines, are available in references 1 to L 
Inasmuch as reliance on such published impingement data is often 
necessary in the design of protective systems, an extension of the use 
of the available two-dimensional (nonswept) impingement data to include 
swept wings is desirable. 

The analysis conducted at the NACA Lewis laboratory to extend the 
use of the two-dimensional impingement data to swept wings is presented 
herein. 
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SYMBOLS 

The following symbols are used in this report: 
a droplet radius 

Cp drag coefficient for spherical droplets in air, dimensionless 

d droplet diameter 

Em collection efficiency, dimensionless 

K inertia parameter, dimensionless 

L cylinder radius, or airfoil chord length (measured normal to lead- 
ing edge) 

Re Reynolds number with respect to droplet, dimensionless 
Reo free-stream Reynolds number, dimensionless 

S distance on surface of airfoil measured from leading-edge chord 
point, ratio to chord length 

T airfoil thickness, ratio to chord length 

t time 

U free-stream velocity 

u dim ensionless local air velocity, ratio of local air velocity to 

either free-stream velocity or projection of free-stream velocity 
in normal plane 

v dimensionless local droplet velocity, ratio of actual droplet 

velocity to either free-stream velocity or projection of free- 
stream velocity in normal plane 

W m rate of total impingement of water per unit span length 

Wp rate of local impingement of water in droplet form 

w liquid-water content 
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x,y,z dimensionless coordinates, ratio of coordinate to either chord 
length or cylinder radius 

a angle of attack 

. P local impingement efficiency, dimensionless 

Y angle of sweep or yaw 

P viscosity of air 

P a density of air 

P w density of water ' 

T time, dimensionless 

scale parameter, dimensionless 
t* stream function 

Subscripts: 

Z lower 

max maximum 

n referred to normal plane 

s referred to free-stream plane 

u upper 

x component in x-direction 

y component in y-direction 

z component in z -direction 

ANALYSIS 
Flow Field 

A knowledge of the flow field about a swept wing is required in 
the droplet -impingement calculations, because the motion of a droplet 
is governed by the momentum of the droplet and by the drag forces 
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imposed on the droplet by the relative motion between the droplet and 
the air which is moving along the streamlines around the wing. The 
droplet momentum tends to keep the droplet moving in a straight path, 
while the drag forces tend to force the droplet to follow the stream- 
lines. 

Because of its geometric and aerodynamic simplicity, the problem 
of impingement on a yawed circular cylinder will be considered first. 

The method developed for the cylinder will apply equally well to a 
swept wing, and this application will be discussed in detail in a sub- 
sequent section. 

The incompressible nonviscous velocity field around a yawed circular 
cylinder of infinite span can be divided into components in the x-, y-, 
and z-directions, as shown in figure 1. The right-rectangular coordinate 
system is established so that the center line of the cylinder coincides 
with the z-axis. The angle of yaw y is measured with respect to the 
negative x-axis. The flow in the xy-plane (normal plane) is the same 
as the flow about a nonyawed circular cylinder which has a free-stream 
velocity equal to the component of the free-stream velocity in the normal 
plane Us cos y. The velocity component in the z-direction is a constant 
throughout the flow field, and the magnitude of the component is deter- 
mined by the projection of the free-stream velocity on the z-axis 
U s sin y. 

The effect of yaw on the streamlines and local air velocity vectors 
can be illustrated by an example in which the flow field about a yawed 
circular cylinder is presented. For this example, the yawed cylinder is 
assumed oriented in an x,y, z- coordinate system as shown in figure 1. 

The angle of yaw y with respect to the free-stream direction is 30°. 

The projection of two streamlines, selected for illustration, upon the 
xy- and xz-planes is shown in figure 2. (The method of calculation is 
given in appendix A.) Also shown are the projections of velocity vec- 
tors upon the xz-plane at various points in the sheet (surface perpen- 
dicular to the xy-plane containing the streamline) which contains each 
streamline. Considerable deviation of the streamlines from the projec- 
tion of the free-stream direction in the xz-plane is indicated in the 
figure. The curvature of the streamlines in the xz-plane is due to the 
fact that, whereas the u x component changes in magnitude with x and 
y, u z is a constant independent of x or y. The z velocity com- 
ponent u z is a constant throughout the flow field, because in a 
frictionless fluid it is equal only to the projection along the z-axis 
of the free-stream velocity U s and is not influenced by the perturba- 
tion velocity of the cylinder in the xy-plane. The maximum curvature 
of the streamline projections in the xz-plane takes place in those 
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streamlines located near a surface defined by the stagnation stream- 
lines. This maximum curvature of the stagnation streamlines is to be 
expected, because along a stagnation line where u x becomes small 
(close to the body), u z is relatively predominant and the displacement 
is primarily in the z-direction. 


Method of Calculating Droplet Trajectories 

Inasmuch as the projection of the path of an air particle in the 
xz-plane does not remain in the free-stream direction but rather is a 
curved path, the projected path of the droplets in the xz-plane will 
also be curved. Because of inertia, the droplets will be deflected 
less from the free-stream direction than will the air. The dimension- 
less equations that describe the motion of droplets along a trajectory 
in a three-dimensional flow field ’can be written in terms of the x-, 
y-, and z-components as follows (refs. 5 and 6): 


where 


dv 

<rr 


X 


dv. 


dT 


y _ 


CpRe 

"24” 

Cj^Re 

~24“ 


T = 


K ( u x “ v x) 
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( 2 ) 


The Reynolds number Re with respect to the droplet is obtained 
from the following relation: 


= t Ux " Vx)2 + (u y ■ v y) 2 + ( u z - v z) : 


( 3 ) 


where 
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The coefficient of drag Cp is a function of Re. The values of 
Cd corresponding to various values of Re are obtained from experi- 
mental drag data (ref . 7 ) . 

The solution of equation (l) to obtain three-dimensional droplet 
trajectories about a yawed cylinder of infinite span is not very much 
more difficult than the two-dimensional solution for the nonyawed 
cylinder, because the air velocity components u x and uy are inde- 
pendent of z and because the droplet and air velocity components in 
the z-direction are equal. Because v z = u z , there is no z -component 
of drag, and equation (lc) reduces to: 


dv z 

dT 


= 0 


and 


z = u z T + constant 

Furthermore, equation (3) reduces to 


(*%;) = (Ux ’ Vx)2 + (Uy " Vy)2 


( 5 ) 

( 6 ) 


Because of the properties of the flow field discussed in the pre- 
ceding paragraph, the projection of the droplet trajectories on the 
xy-plane can be obtained by the simultaneous solution of equations (la) 
and (lb) with the aid of a differential analyzer. The manner of cal- 
culating is the same as the two-dimensional method presented in refer- 
ence 6, if the component of the free-stream velocity in the normal plane 
is used throughout the computations instead of the free-stream veloc- 
ity U s . In order to perform this calculation, both sides of equa- 
tions (l) and (6) are divided by cos^. Thus, all velocities are used 
in the form of a ratio to the normal component of the free-stream veloc- 
ity) T becomes T n = tUn/L) and equations (2), (4), and (6) become 
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respectively. The di m ensionless equations of the x- and y-components of 
motion of the droplet along a trajectory are now in a form identical to 
that of the equations for a nonyawed cylinder. In addition, the x- and 
y-components of the velocity field about a yawed cylinder will be in a 
dimensionless form identical to that of the x- and y-components for a 
nonyawed cylinder (appendix A) . 

The projection of the droplet trajectories upon the xz-plane is 
obtained by the simultaneous solution of equations (la) and (5), which 
is accomplished if the relation between x and T n is known. This 
relation between x and T n can be obtained from a differential 
analyzer as it solves the xy-component of the trajectories. The value 
of z corresponding to the various values of x can be obtained from 
equation (5), where u z ^ n = v z ^ n = tan r (appendix A). 


Sample Droplet Trajectories About Yawed Cylinder 

Three-dimensional droplet trajectories were calculated for a yawed 
cylinder for conditions where the inertia parameter K n is equal to 
2.78 and the free-stream Reynolds number in the normal plane Reo n is 
equal to 100. The projections upon the xy- and xz -planes of typical 
droplet trajectories are shown in figure 3. The starting ordinate 
(at x = -co) in the xy-plane for each trajectory is given by the value 
of y 0 in the figure. 


Swept Wings 

The calculation of droplet trajectories about a swept wing can be 
accomplished by the method used for the yawed cylinder in the preceding 
sections, if the flow field about a swept wing is approximated (except 
near the fuselage and wing tips) by the incompressible non viscous flow 
about an infinite-span wing of zero taper at an angle of yaw equal 
to the sweep angle. This approximation is valid for wings of high 
aspect ratio and small taper. 

The nonviscous flow field about a yawed wing of infinite span can 
be .obtained in the same manner as that for the yawed cylinder by deter- 
mining the components in the normal (xy) plane (fig. 4) and the com- 
ponent parallel to the leading edge (z-direction) . The flow field in 
the normal plane is determined from two-dimensional theory by using the 
airfoil section in the normal plane and the component of the free- 
stream velocity in the normal plane U s cos r. As before, the 
z-component is constant throughout the flow field, being simply the 
projection of the free-stream velocity upon the z-axis U s sin r. 
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When the swept wing is at a free-stream angle of attack a s 
(fig. 4), the effective angle of attack in the normal plane, to a close 
approximation, is <% = cig sec T and the components of the free-stream 
velocity are U s sin r a and U s cos r a , where r a is the angle between 
U s and U n . For small angles of attack, r a is for all practical pur- 
poses equal to r and therefore r can be used when calculating the 
components of the free-stream velocity; however, the corrected angle of 
attack a s sec T must be used in obtaining the flow field in the normal 
plane . 

Air streamlines about a 65i-212 airfoil at 45° yaw and an angle of 
attack in the normal plane of 4° (a s = 2.8°) are shown in figure 5. 

The projection of three streamlines on the xy-plane from x = -1 to 0.1 
is shown in figure 5(a), and the projection on the xz-plane is shown in 
figure 5(b). The method used to obtain the streamlines is given in 
appendix B. As for the cylinder, the projections of the streamlines in 
the xz-plane deviate from the projection of the free-stream direction 
because of changes in u x . The deviation from the projection of the 
free-stream direction is considerably less for the airfoil than for a 
cylinder with the same chord length, as would be expected from compari- 
son of the velocity field of a thin body with that of a blunt body. 


Sample Droplet Trajectories About Yawed Wing 

The projections upon the xy- and xz-planes of typical droplet 
trajectories about a 65p-212 airfoil with a s = 2.8° (tin = 4°) and at 
a yaw angle of 45° are shown in figure 5 for values of K n = 0.1 

and Re 0 ^ n = 128. As in the case of the air streamlines, there is less 
curvature of the xz-projection of the trajectory for the airfoil than 
for the relatively blunt cylinder. For example, up to about 1 chord 
length (x = -1, fig. 5(b)) ahead of the wing the droplets are essen- 
tially undeviated from the free-stream plane; whereas, for the cylinder 
they are undeviated only up to about 2 diameters (x = -5, fig. 3) ahead 
of the cylinder. 


RESULTS AMD DISCUSSION 

It was shown in the ANALYSIS that the equations of the x- and 
y-components of motion of the droplet trajectories and of the air- 
velocity field around a swept wing can be put in the same dimensionless 
form as the x- and y-components of motion for the two-dimensional 
(nonswept) case. Because of the identical form of the dimensionless 
equations, it is possible to obtain the xy-projection of the droplet 
trajectories around a swept wing from two-dimensional trajectory data 
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by using the component of the free-stream velocity in the normal plane 
and the angle of attack in the normal plane corresponding to the given 
angle of sweep T when evaluating the various dimensionless parameters. 
The projection of the trajectories in the xz-plane cannot be obtained 
from two-dimensional data without recalculation with a differential 
analyzer. A knowledge of the projection of the trajectories in the 
xz-plane is not usually needed in order to obtain an appraisal of the 
impingement of droplets on a swept wing, because the motion of the drop- 
lets in the xz-plane does not influence the chordwise droplet- 
impingement distribution calculated from trajectories in the xy-plane. 

The spanwise displacement of the droplets does cause a small spanwise 
shift in the point of impingement. For droplets of uniform size, the 
xy-impingement pattern is thus displaced spanwise a small amount; but 
as the trajectories in the xy-plane are not dependent on z, the impinge- 
ment distribution on an airfoil section in the normal plane is the same 
at all spanwise (z) positions along the wing. If the cloud consists of 
droplets having a distribution of droplet sizes, each size will be 
shifted spanwise a slightly different amount; but the net chordwise 
impingement in the normal plane at each spanwise station will be the 
same as the net impingement calculated from the xy- component of the 
trajectories when the given droplet-size distribution is used. 

Example of application of two-dimensional data to swept wing . - 
The determination of droplet impingement on a swept wing from two- 
dimensional data calculated for nonswept wings, such as presented in 
references 1 to 4, is best illustrated with an example. 

For the example, a wing with a 65^-212 airfoil section in the normal 
plane is swept to an angle of 45°. The chord length in the normal plane 
is 8 ^ feet, and the free-stream angle of attack is 2.8° (a^ = 4°). The 

swept wing is assumed to be moving in the free-stream direction at 
525 miles per hour through a cloud composed of droplets all of which are 
20 microns in diameter; the liquid-water content of the cloud is assumed 
to be 0.2 gram per cubic meter; and the free-air temperature is -11° F 
at a pressure altitude of 20,000 feet. 

The extent of impingement on the airfoil, the rate of collection of 
total water, and the rate of impingement per unit area of airfoil surface 
can be found for this swept wing from two-dimensional data for an MCA 
65^-212 nonswept wing presented in reference 2, provided the values of 

the dimensionless parameters required for entering the data figures of 
reference 2 are calculated for the normal plane. The required dimension- 
less parameters K^, i|r n , and ReQ^ n * >or the swept wing are 
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U n = U s cos T = 371 mph 

d 2 U 

K n = 1.704X10" 12 — =— = 0.09 
u |oL 


dp U 

Re n „ = 4.813X10 -6 - = 131 

[X 

and 


r n 


Re 0,n 

*n 


= 1455 


(Factors required for converting from self-consistent units required in 
eqs. (l) to (9) to units most commonly used in aeronautics are given in 
ref. 2.) Typical droplet trajectories which approximately correspond 
to these values of K n and Reg^ n are given in figure 5. From the 
method of reference 2, the farthest point of impingement on the upper 
surface ^ of the. airfoil obtained from figure 6(a) (\(r = \|r n and 

Reg = Reg n ) is 0.025 chord length, and the farthest point of impinge- 
ment on the lower surface (fig. 6(b)) is 0.14. The collection 

efficiency corresponding to these values of Reg n and \|r n 

is 0.24 (fig. 7). The rate of total water in droplet form impinging 
on the swept wing in pounds per hour per unit span (parallel to the 
leading edge) is 


W m = 0 . 329E m TnJ s w cos T = 5.86 lb/(hr)(ft span) 

The term cos T is introduced into the equation to account for the 
increased spanwise length of wing sweeping out a given volume of cloud 
in comparison with an identical nonswept wing moving with the same free- 
stream velocity. That is, the water collection rate per unit span is 
reduced by cos T for the swept wing, which is consistent with the 
mathematical artifice of considering the airfoil to be moving forward 
in the normal direction with a velocity U s cos Y". 

The maximum local impingement efficiency factor obtained 

from figure 8 is 0.55, which occurs at 0.008 chord length along the 
lower surface. The location of the region of maximum local impingement 
is of importance when determining the most effective location of the 
parting strip on a cyclical thermal de-icing system. The maximum rate 
of local impingement is 


W P,max = °* 329u s w Pmax cos “n cos T ~ l 3 * 39 lb/(hr)(sq ft) 
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Similar calculations for an identical nonyawed wing moving at the 
same free-stream velocity (525 mph) in an identical environment 
(d, 20 microns; w, 0.2 g/cu m; temperature, -11° F at 20,000 ft pressure 
altitude) give the farthest point of impingement on the upper surface as 
0.028 chord length and the farthest point of impingement on the lower 
surface as 0.17 chord length, the collection efficiency Em equal to 
0.28, the total water catch as 9.67 pounds per hour per foot span, the 
maximum local impingement efficiency (3 mav equal to 0.58, and the maxi- 
mum rate of local impingement Wp jiaax equal to 19.98 pounds per hour 
per square foot. (For purposes of comparison, the assumption has been 
made in this paragraph that the impingement data of reference 2 are 
valid beyond the critical Mach number of the nonswept wing, which may 
not be true. ) 

The method for obtaining the impingement on a swept wing was derived 
for a wing of zero taper. The calculations for a wing with small taper 
can be made from the impingement data in the same manner as suggested 
in reference 2 for nonswept wings with small taper. 

The impingement on a yawed cylinder can be obtained from two- 
dimensional data by a procedure similar to that illustrated for the 
swept wing. The pertinent dimensionless parameters, conversion factors, 
and practical units for the circular cylinder are given in reference 6. 

Transformation of airfoil section from free-stream plane to normal 
plane . - The method discussed for obtaining impingement on a swept wing 
from two-dimensional data requires a knowledge of the airfoil section 
in the normal plane. If a swept-wing layout has its design airfoil 
section in the free-stream plane rather than in the normal plane, then 
it is necessary to determine the airfoil section in the normal plane. 

This transformation can be made by multiplying the airfoil ordinates 
by the secant of the sweep angle. Thus, the airfoil section found for 
the normal plane will have a thickness ratio which is larger by the 
factor sec X. For example, at 40° sweep angle a wing with a 6-percent- 
thick airfoil section in the free-stream plane will have approximately 
an 8-percent-thick (7.83) airfoil section in the normal plane. In 
addition, for most airfoil types the section in the normal plane will 
be, to a satisfactory approximation, a member of the same airfoil 
family as the airfoil section in the free-stream plane. Noncambered 
airfoils generally transform in this manner so as to correspond more 
closely to a member in the same family than do those with appreciable 
camber. 

Value of xz-projection of droplet trajectories . - Although the 

projections of the droplet trajectori e s upon the, xz-p.lane. are., not .needed 

for determining droplet impingement upon a swept wing, a knowledge of 
the xz-projection is of value when considering the impingement of drop- 
lets on small wing appendages. This knowledge is particularly useful 
when considering the impingement on appendages located near the leading 
edge of a swept wing, such as boom-mounted pres sure -sensing elements, wing 
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fences, and small leading-edge inletB. Although these devices will 
change the local flow field to some extent, the impingement zone would 
generally be expected to extend farther back on the inboard side of 
protuberances and on the outboard side of cavities than on the opposite 
or "shadowed" side of the appendages attached to a sweptback wing. 

For appendages attached to a sweptforward wing the same thing would 
occur, except that "outboard" and "inboard" should be interchanged in 
the preceding sentence. 

Application to supersonic velocities . - The preceding sections have 
specifically treated only the extension of droplet -impingement data for 
nonswept wings to swept wings moving at subsonic velocities in an 
incompressible fluid. The same general approach can be employed to 
obtain impingement data for swept wings of high aspect ratio moving at 
supersonic speeds from nonswept-wing impingement data calculated for 
supersonic velocities in a compressible fluid. This extension is 
possible because, irrespective of the magnitude of the velocity, the 
component of translation of any cylindrical body in the direction of 
its long axis has no effect on the motion of a frictionless fluid. The 
flow around a wing of constant section moving through a fluid at rest 
is determined by the normal components of velocity of its solid bound- 
aries, and these components in turn are completely specified by the 
component of motion in planes perpendicular to the axis Us cos T 
(ref. 8). Thus, the impingement on a swept wing of high aspect ratio 
moving at supersonic velocity can be determined from two-dimensional 
(nonswept-wing) data calculated for a wing with the same airfoil sec- 
tion as the normal airfoil section of the swept wing and moving through 
an identical environment at a Mach number equal to the Mach number found 
for the normal plane of the swept wing. 


CONCLUSIONS 

An analysis of the general effect of sweep on the droplet trajec- 
tories about a swept airfoil arid a method of calculating the droplet 
trajectories have been presented. The extent of impingement of cloud 
droplets on the airfoil surface, the rate of collection of total water 
in droplet form, and the local rate of impingement per unit area of the 
airfoil surface can be found for a swept wing from two-dimensional data 
for a nonswept wing with the same airfoil section as the section in the 
normal plane of the swept wing. The values of the dimensionless param- 
eters required in order to interpret the two-dimensional data are cal- 
culated with respect to the normal plane of the swept wing rather than 
with respect to the free-stream plane as is the usual procedure for 
nonswept wings. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 12, 1953 
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APPENDIX A 


CALCULATION OF STREAMLINES AND VELOCITY VECTORS ABOUT 
CIRCULAR CYLINDER AT YAW ANGLE T 

A cylinder oriented with the axis along the z-axis of an 
x,y, z-cartesian coordinate system (fig. l) is considered for the three- 
dimensional derivation of the velocity field around a yawed cylinder. 

The stream function for nonviscous incompressible flow in the x'y' -plane 
using polar coordinates r* and 0 and letting i|r* = 0 at the surface 
of the cylinder is 

** = U n (r* - ^ sin 0 (Al) 

where L is the radius of the cylinder, and the terms with the prime 
superscript have the dimensions of length. Or, in dimensionless form 

iy; - ( r - f) sin e < A2 > 

where (r = r'/L) and ijr*/u n L are dimensionless. For various i|r*/u n L 
values, the x- and y-coordinates (x and y are dimensionless), which 
can be obtained from equation (A2), give the projection of the stream- 
line upon the xy-plane. 

The, projection of the streamline in the xz-plane can be obtained 
as follows: The streamlines in the xy-plane are divided into small 

line segments. The x-component of velocity corresponding to the 
i th segment is determined from the average of the velocities determined 
by the x and y values of the end points of the i th segment. The 
time required for an air particle to travel the ith segment of a trajec- 
tory is 


(AT n )i 


**i 

( u x, n)i 


The displacement in the z-direction while the air particle travels 
a distance Ax^ in the x-direction is 


^i = ( AT n)i K,n)i 
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From any arbitrary point ahead of the cylinder, for each Xj_ 

i i 

z i = 2, (^ T n)i ( u z,n)i = u z,n 2 ^ T n)i 
0 0 

since u z ^ n is constant. 

The dimensionless velocity components of the flow field about a 
circular cylinder at yaw angle X are 


u x 

Ux > n “ U s cos r 


x 2 . y 2 
(x 2 + y 2 ) 2 


- ^ 

l Vj n U s cos X 


2xy 

(x 2 + y 2 ) 2 


u z U s sin r 

- u s cos r = u s cos r = tan r 

where the velocity components with the prime superscript have the usual 
dimensions associated with velocity. 

Projections upon the xy-plane of streamlines about a 30°-yawed 
cylinder for \|r*/u n L = 0.19 and i|r*/u n L — 0 (infinitesimally greater 
than zero) are shown in figure 2(a). The corresponding projections 
upon the xz-plane along with the projections in the xz-plane of the 
velocity vectors in the sheets containing the selected streamlines are 
shown in figure 2 (b ) . 
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APPENDIX B 


CALCULATION OF STREAMLINES AND VELOCITY VECTORS ABOUT WING WITH 
65 1 -212 AIRFOIL AT YAW ANGLE T AND ANGLE OF ATTACK a g 

For a wing in an x,y,z- coordinate system with its leading edge 
along the z-axis and the free- stream direction at an angle Y with 
the negative x-axis (fig. 4), the equations of the velocity components 
of the air flow about the wing can be written (by using the vortex 
method of ref. 2 and s ummi ng the contributions of 300 vortex elements 
distributed on both the upper and lower surfaces of the airfoil.) 


u x = U s cos 


300 

Z 

i=0 


u s (y-T)) 


2jtr" 


AS + U c 


COS r cos 0^ 


u' = U,, cos 
y s 


300 

Z 

i=0 


u s(&- x ) 


AS + U 


2itr" 


s 


cos Y sin 


u* = U sin T 
z s 

where \ and t) are dimensionless coordinate points (ratio to chord 
length) on the airfoil, r is the distance (ratio to chord length) from 
an element of the vortex sheet to a point in the flow field, S is the 
distance (ratio to chord length) from chord line measured along the 
surface of the airfoil, the subscript S refers to the airfoil surface, 
and the prime superscript indicates that a term has the usual dimen- 
sions. 

Dividing both sides by the normal velocity component yields 


u' 


u 


x,n 


U. s cos T 


u * 


u n - 


300 


i=0 


Ug(y-ii) 

2itr 2 


AS + cos 
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300 

u* u» \ — ' 

y y \ 

u y> n s u s cos r s Un7 , 

i=0 


u s (e-x) 

2 — AS + sin (% 


2itr £ 


u z u z 

= rT- = tan T 


z > n U s cos T " U n 


The summation terms in the equations were computed with the aid of a 
punched-card electronic calculating machine (ref. 2). 

Inasmuch as a closed-form solution of the equation of the stream- 
lines was not available, the streamlines were obtained by a simultaneous 
solution of the following equations with the aid of the differential 
analyzer: 


x = *0 + / "x,n dT n 

y = y 0 + / Uy,n dT n 

z = z 0 + u z,n T n 

where the subscript 0 refers to the starting coordinate. 

The projections of the streamlines in the xy- and xz -planes for 
a 65^-212 airfoil at a normal angle of attack of 4° (og = 2.8°) and 
45° angle of yaw are shown in figure 5. 
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(b) Projection upon xz-plane. (Arrows represent 
of local air velocities in each sheet 








Figure 4. - Geometry applicable to yawed wing. 
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(a) Projection upon xy-plane. (Ordinate scale expanded x4.) 

Figure 5. - Typical air streamlines and water-droplet trajectories in the vicinity of 65^-212 airfoil at yaw angle of 45°. Angle of attack 
a nj 4°; inertia parameter K n , 0.1; free-stream Reynolds number Reo,n> 128. 
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(b) Projection upon xz-plane. 

Figure 5. - Concluded. Typical air streamlines and water-droplet trajectories in 
the vicinity of 65 1 -212 airfoil at yaw angle of 45°. Angle of attack 4°; 

Inertia parameter K n , 0.1; free-stream Reynolds number Re 0 n , 128. 
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Figure 6. - Farthest point of droplet iu^ingement on airfoil surface. NACA 65^-212 airfoil; angle of attack, 
(Reproduced from fig. 6, ref. 2.) 
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Figure 6. - Concluded. Farthest point of droplet impingement on airfoil surface. NACA 65^-212 airfoil; angle of attack 
(Reproduced from fig. 6, ref. 2.) 
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Distance on airfoil surface, S, ratio to chord length Upper 

(a) Pree-atream Reynolds number with roBpect to droplet, 16, 
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il Impingement efficiency. NACA 65j-212 airfoil; angle of attack, 4 °. (Reproduced from fig. 10, 
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(c) Pree-atream Reynolds number with respect to droplet, 1024. 

Figure 8, - Concluded. Local impingement efficiency. NACA 85^212 airfoil; angle of attack, 4°. (Reproduced from fig. 10, ref. 



